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Abstract: The "Universality Theorem" for gravity shows that f{R) theories (in their metric-affine for- 
mulation) in vacuum are dynamically equivalent to vacuum Einstein equations with suitable cosmological 
constants. This holds true for a generic (i.e. except sporadic degenerate cases) analytic function f{R) and 
standard gravity without cosmological constant is reproduced if / is the identity fimction (i.e. f{R) = R). 
The theorem is here extended introducing in dimension 4 a 1-parameter family of invariants '-?R inspired by 
the Barbero-Immirzi formulation of GR (which in the Euclidean sector includes also selfdual formulation) . 
It will be proven that f(^R) theories so defined arc dynamically equivalent to the corresponding metric- 
affine f(R) theory. In particular for the function f(R) = R the standard equivalence between GR and 
Hoist Lagrangian is obtained. 



1. Introduction 

It is well-known that for almost any analytic function /(i?) the metric-aiiine theory with 
Lagrangian Lf{g,j^r) = \/g f {R{9 1 j^^)) is dynamically equivalent to standard GR with a 
suitably quantized cosmological constant (encoded by the function /); see [1] and [2]. Hereafter 
j^r refers to the fact that the Lagrangian depends on the connection V and its first derivatives. 

The original universality theorem was established in vacuum in [1]. In fact, matter coupling 
produces in /(i?) theories new effects with respect to Einstein equations, since a new metric 
conformal to the original one can be defined (see [3]). The theory in these new variables is quite 
similar to standard GR though with the addition of an effective energy momentum tensor; see 
[2]. Such new effects have been recently investigated aiming to find a specific f{R) theory able 
to model dark energy and dark matter phenomenology; see e.g. [4], [5], [6], [7]. 

The metric-affine formulation has been recently criticized. The most serious criticism is based 
on a theorem which shows that, under specific state equation hypotheses, internal and external 
solutions for a (stationary, spherically symmetric) polytropic star do not match on the boundary 
and produce singularities on the star surface; see [8]. However, as it often happens when precise 
no-go theorems are formulated, it has been later argued that the hypotheses of this theorem 
are in fact physically unreasonable; see [9]. In this particular case it has been shown that (at 
least for specific examples of /) the singularity depends on matter densities far lower than the 
one that is reasonably expected. We shall not discuss further these aspects here. 

Another criticism to /(i?) theories (this time in purely- metric formulation; see [2]) is based 
on fixing the variation of first derivatives of field on the boundary; see [10]. This criticism 
seems to be based on a physically and also mathematically questionable method discussed 
[11], that is meaningful only in standard GR; in standard purely-metric GR, one can in fact 
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subtract a suitable boundary term from the action so that fixing first derivatives of the metric 
on the boundary (as it is prescribed by Calculus of Variations; [12]) is not necessary any longer 
provided one accepts to modify the Hilbert Lagrangian by adding suitable ad-hoc boundary 
counterterms. The same unconventional procedure, however, cannot be carried over for a generic 
f{R). We argue that one cannot consider this a problem; see [13] for a detailed criticism of 
the method. In standard (higher order) variational calculus all derivatives of variations of fields 
are fixed up to one order less with respect to the effective order of the theory. The fact that 
first derivatives do not need to be fixed in GR, which is a second order gravitational theory 
when formulated in terms of the Hilbert Lagrangian, is just a consequence of the well-known 
fact that GR is in fact degenerate and can be formulated even covariantly as an equivalent first 
order theory; see [14]. Moreover, not fixing higher order variations causes unreasonable results 
in many cases, but wc shall not discuss these aspects here, either; sec [13]. 

We shall here establish the dynamical equivalence between f{R) theory d la Palatini and its 
Barbero-Immirzi formulation. The special case for f{R) = R is already well-known; see [15], 
[16], [17], [18], [19]. The extension to a generic f{R) is new, to the best of our knowledge. 

2. Notation and Hoist Formulation of GR 

Hoist formulation of GR is the classical basis of LQG formulation in terms of the Barbero- 
Immirzi connection. 

Let us consider a 4 dimensional (paracompact, connected, orientable) manifold M which allows 
metrics in signature rj = (3, 1) and global spin structures (i.e. with zero 1st and 2nd Stiefel- 
Whitney classes). Let P be a principal bundle over M with SO{r]) as structure group. Notation 
follows [20]. 

Let us denote by J^P the first jet prolongation of P; there is a standard right action of S0(?7) 
on J^P, namely the prolongation of the canonical right action Rg : P ^ P on P; see [20]. The 
bundle C(P) = J^P/SO{ri) admits a global family of local coordinates (x^jV^) and (global) 
sections are by definition (global) SO(77)-connections on P. 

Let A : (GL(4) x SO{t])) x GL(4) ^ GL(4) : ( J, ^, e) J • e • be the natural action of the 
group GL(4)xSO(r7) on the manifold GL(4). The associated bundle e(P) = (L(M) xP) XaGL(4) 
has coordinates {x^,e^) and (global) sections, that always exist under our hypotheses on M, 
are by definition (global) frames. 

Here frame is meant in a non-standard sense. A global frame for us is not a global seetion of the 
frame bundle L{M), that may not exist, but rather a global section of e(P). A global frame is 
here a family of local sections of L(M) defined on an open covering of M and such that transition 
functions arc valued in S0())) rather than in {II} as required for global sections of L(M). 
This is done since generic manifolds do not allow global sections of L{M); the topological condition 
for this is physically too strong in general, while e(P) always allows global frames in the above 
hypotheses; sec [20]. 

Moreover, the standard covariant derivative of sections of e(P) is the covariant derivative introduced 
ad hoc and used in literature for frames, which is not the standard covariant derivative for sections 
of L(M); see [20]. 
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Vielbein ej^ are frame inverses; vielbein induce a metric on M 

gi,u = elr]abel (2.1) 

where jy"'' is the standard diagonal matrix of signature jy = (3, 1). Hereafter Greek indices are 
moved up and down by the metric g^i, while Latin indices are moved by rjab- 

Let us now define the curvature 2-form 



j^ab ^ Ij^ab^^ A dx" (2.2) 

where R"'' is the curvature tensor of r"*, the local representative of an arbitrary global 
section F of C{P). The vielbein form is the 1-form 

e" = el dx" (2.3) 



For the Hoist formulation (see [17], [19], [21]) let us set C = C{P)x Me{P) for the configuration 
bundle; the Lagrangian is: 

LH = i-, {R"" a a e'eabcd - ^R"' A e„ A e,) (2.4) 
where k and 7 are constants. For later covanience this Lagrangian can be written as 

, / (2.5) 

where ds is the standard local volume form induced by coordinates, e is the (module of) the 
determinant of the vielbein e^ and we set /3 = — 5^ and 

I'R = R'\.e''A + m^\.e'''e'"eM (2.6) 

Let us also set R := R'^^^^e'^e^ and i?^ := R'*^„e'^ which are functions of (e, j^F). 

Let us stress that the additional term R""^ fw^''''' e''-" tcdab in ^R is peculiar of the Palatini formu- 
lation; if we assumed the connection to be metric from the beginning, the curvature would be 
the Riemann tensor of a metric and the additional term would vanish identically due to the 
symmetry properties of Riemann tensors, in particular the first Bianchi identity R"\p^v\ = 0. 
One obtains the same result also using no frames, but allowing torsion; it is sufficient to use 
the correction e^'^''°'Riij,p„{j^r) which does not vanish identically if torsion is allowed. Here we 
have chosen to use frames in view of possible applications to LQG or possible applications with 
spinor couplings. 

We refer to [17] or [21] for the equivalence between the Lagrangian (2.5) and standard GR; 
here it will follow from the universality result proved below, in the particular case f{^R) = ^R. 

Before proceeding to consider these further extended gravitational theories let us first briefly 
review the standard results about mctric-affine f{R) theories. The metric-affine formulation of 
f{R) models is described by the following Lagrangian: 

Lf{9,j'T,j'<f>) = ^f{R) + Cmai{9,3'<t>) (2.7) 
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where F is now a (torsionlciss) linear connection and is a set of matter fields; the matter 
Lagrangian does not depend on F (though it could depend on the Levi-Civita connection of the 
metric g). 

The variation of the Lagrangian Lf with respect to the metric g and the connection T is: 

SLf =^g {-~\fg,Jg"" + f {5g'-'''R,u + g'"'5R,,)) - ^gTi^Jg"" = 

=V9 {-\f9,n.5g>"' + /' {5g""R,, + 2g"'^^ y^Sul^)) - ^T„Jg>-'^ = (2.8) 

where f and f = -tb arc understood to be evaluated at R and we set u^„, = T^,,, — r"> 

•' dR III' (p v)a 

and — y^T^y = One should also consider variations with respect to matter fields, which 

account for matter field equations; we are not interested here in that. 

As usual the last boundary term in the variation vanishes because of boundary conditions 
(5u^^ = since SV^^ = 0) and field equations are 

{/ Rlpv) — If duv = 
Va iVgfgn = 

By tracing the first equation with g, one obtains the master equation f'{R)R — 2/(i?) = T, 
being T = g^^T^j^i,; except in degenerate cases this can be solved for R{T). 

Let us now set a = sign(/') and define a new metric gij,^ = \f'\gij,v; accordingly one has 
a^/gg'^" = ^/gf'g'^"- This can be used in the second field equation which implies that F^j, are 
nothing but the Christoffel symbols of the metric g. 

We can replace this information back into the first field equation and obtain 

f R^.-\f9^. = T^. (2.10) 

where R^i, is now the Ricci tensor of the metric g. Let us stress that /' and / are still evaluated 
along the Ricci scalar R of the original metric g. 

This equation can be recasted as an Einstein equation for the conformal metric g 

Gfii/ := Ri^v — ^RQi^i/ = jr (Tuv — \ {{f)^^~ /) 3/ti') = jrT^i/ + (2-11) 

with an additional eflfective source T^J = {{i')'^R^ 9p.v- The conservation laws for 

matter follow from Bianchi identities of g and read as follows 

{^jrTl"' + {T^S)^^'^ = (2.12) 
where V is the covariant derivative with respect to g. 

Lot us stress that the offectivc energy-momentum tensor T^J , being effective, does not need to 
obey other separated physical energy conditions which are usually required for matter energy- 
momentun tensors. In this sense one can hope to choose / so that the effective energy-momentum 
tensor mimick dark energy/matter without needing to introduce extra exotic matter fields. Of 
course, the freedom in choosing matter fields to model the dark side of the universe is transformed 
into the freedom in choosing the function /. 

The standard universality result has been proved in the mctric-affiuc formulation (2.7); it holds in 
vacuum (i.e. Tf^y = and, more generally, when the trace T vanishes); the master equation reads 
thence as 

f'R-2f = (2.13) 
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which, except few particuhir cases and the degenerate case f = B? that makes it an identity, has 
isolated zeros; by choosing po to be one of the simple zeroes then one must have R = pQ on-shell 
and the field equation can be recasted as 

/' - 3-R5|Ui/) =0 ^ R^i, - \kg^v = --^p]^PQ9nu = Ag^i^ (2-14) 

which are in fact Einstein equations with cosmological constant A = — 4|j7yPo- Multiple zeroes are 
discussed in [1]. 



3. Universality Theorem 

Let us now consider the Lagrangian 

L = e/(^i?)+L™t(e,jV) (3-1) 

where e = det(e°) denotes the determinant of the frame, / is an analytic function such that the 

so-called master equation ^Rf'{^R) — 2f{^R) = T can be solved for ^R{T) as in Section 2. Here 
(/) is a set of matter fields; the important issue is the hypothesis that the matter Lagrangian is 
independent of F^'' (while it could depend on the spin connection oj"^''{e) induced by the frame). 

The variation of the matter Lagrangian defines the matter energy-momentum tensor 

-2eT; = ^ (3.2) 

Let us start computing the variation of the Lagrangian (3.1) with respect to the frame and 
the connection: 

' ' 



SL=e[- fe;de>i + /' (2 V^^rf ej^e^ + 2i^^^e^ 

=2e {f'R-, - If el + pR'^A^od"' - T^) ^e^+ 

- 2V, (e/'e^:^) (5lX} + Pe'^'ab) + 2V, [ef'e^.e^, (dlX] + fi^'-t,) 



(3.3) 



As usual the boundary term vanishes because of the boundary conditions SV^^ = and field 
equations are 

f'R^f, - \ fel + PR'^'i^vGh^cd"^'- = T'^ 

V,(e/'e;^'e;^)(<5[„<]+/3e^:l,) =0 ^^"^^ 

Let us set $1^'^ := 5yX] + P^"''-ah] being it skew in both the upper and lower pair it defines an 
endomorphisms $ : — ^ in the space of skew 2-tensors. It can be proved (e.g. by explicit 
computation by Maple) that it is invertible. One can easily check that its inverse is 

($-i):J = (l + 4/32)(5f,^^-/3e«^ (3.5) 



r 
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Here we are discussing Lorentzian sector. In Euclidean signature the special cases /3 = must be 
dealt with separately. They correspond to (anti)selfdual cases in which the map <&^^ is degenerate. 
Accordingly, the analysis needs extra care though it leads to similar results. 

L 

The second field equation is then 

V.(e/'ef,e:^])=0 (3.6) 

Let us now define a new frame 



e 



WK (3.7) 



We shall systematically denote by a tilde the quantities computed in the new frame. 
The second field equation is then 

v. (e-gf.g;^) = (3.8) 

which is the same equation obtained in the standard frame-affine framework; it implies that 
T^^ = Cj"^ coincides with the spin connection induced by the frame ef^. 
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Equation (3.8) can be recasted as: 

V^(ee^fe^6'"'«<i)=0 ^ V;, (eS^e^/^M-) = ^ (e^e^,) = (3.9) 

Let us define fc^^'c = {t'^ — (^^^ Cc ; this is the difi^erence between the dynamical connection T 
and the spin connection induced by the frame e. By construction k"'^ = —k^"^. By considering 
the covariant derivative of the frame we have 

e^V^e« = (<i^e« - f ^^e^ + T\4) = (v^e«) + k%ei = k%4 (3.10) 

Hence equation (3.9) can be recasted as 

<^ k"egS) + k"'j,/g + k"'gjSl = kK.gS'j + k^'fj"g + k^gfSi (3-11) 

=^ Zk^eg "f" '^^^ ge — ^ gb^e ^ po, — 4fc^pa ^ gii ~ ^ 
and substituting back into the original equation we obtain 

Zk"-eg + Sfe^ge = ^ k^eg = -2fe''(ge) ^ ^"[^3] = (3.12) 

Finally, we have the following lemma: 

Lemma: if k^hc = -^bac and k^bc = Kcb t^ien k^bc = 
Proof: Let us simply notice that 

^abc ~ ~^bac ~ ~^bca ~ ^cba ~ ^cab ~ ~^acb ~ ~^abc (3.13) 
from which the thesis follows. ■ 
Hence since fc"*^ = then T'^^ = uf. 
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L 

This piece of information can be used to express the curvature tensor in terms of the Riemann 
tensor of the metric g^^ induced by the frame e^: 

jyab pa Z"- zi 

-O- fJ.1/ — ^ Xiii/y ^a^fi 

R\, = R"\,,el = ^\R\,.~g^f'~e%el = ^{-g^f'Rp.el = ef'r^Rp^el = R%el (3.14) 
R = R^ei^ = ef'R 

(Let us stress that in our notation Greek indices are moved by the metric g and not by g; hence 
R"i,, means R" i, = not g"''^i?fj,j as usual. This is a consequence of having two metrics 

around. One should specify which metric is involved at any step!) 

Before being ready to manipulate the first field equation we have to prove that the extra term 
in equation (3.4) vanishes: 

R'%.e-,ec/.' =VmR%,.g'^elef,eU,X''. = ^AF\R" x^.g'"' e^^e^r ecdebri""' = 

= xAF\eR''x,^e;ea.6p.~g'''f''r = " VlFlie)-' R" x,.e;~9ape^''"' = (3-15) 

= - V\Fm-'~g,ae; R^ox/'^-" = 

and that the scalar ^R coindices with the Ricci scalar R of the metric g induced by the original 
frame eit: 



PR =R + pR^^^e^^e'^'-e^, = R + e^f R''x,.~9^^el^pe'^peUcdahr'r = 
=R + ep~efR"x,,uep,a&rT~9^^ = R- e^(e)- V'^"[AHff«^e'''''^ = ^ 



(3.16) 



The first field equation in (3.4) can be now recasted as follows 

(e{ffg"^Rp,-\f5i)el = T; (3.17) 

- ys;) = T;el =: T% (3.18) 

f'R^, - =■■ T^. (3.19) 

which coincides with the equation in standard metric-affine /(i?)-theories (see equation (2.10)). 
In fact, /' and / arc evaluated at the Ricci scalar R of the original metric g and if the matter 
Lagrangian depends on the frame through its associated metric (as assumed otherwise there is 
no metric aflSne formulation to compare with) one has 

-2eT; = ^ = 2^e- = -l^gT^^e'^'' T,, = T^e^. (3.20) 

We stress that equivalence holds both in vacuum and in the presence of matter. 
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(4.2) 



4. Conservation Laws 

We shall here compute and discuss conservation laws for f{^R) theories, following [20]. We 
shall use the formalism introduced in [22] to which we refer for motivations. The case of 
standard Hoist action has been already discussed in [23] . 

The Lagrangian (3.1) is gauge-natural (sec [20]); any generator S = ^''(9^ + ^""^(Jab of auto- 
morphisms on P is thence a symmetry. Here aab is a right invariant pointwise basis of vertical 
vector fields on P. Accordingly, (in vacuum) the following Noether current is conserved 

£ = e (2/'eJfe^$^i:Hr«^ - ^f) ds, (4.1) 

where rfs^u is the local basis of 3-forms on M induced by coordinates and £e^1^ = ^'^R'^^v/j, + 
"^^^^^^ denotes the Lie derivative and we set ^^^^ = + ^i^Vf for the vertical part of the 

symmetry generator. Hence we have 

=e {e (2/'eJfe^i?«'A. + 2/3/'e^e^e<:<'„i,i?"\. - S^f) + 2/'eJfe^$^V^er')) = 
= {2ee {f'R'x - kfel + PfeAR'Kx) - 2V , (e/'ejfe^) ^^^(t)) ds,+ 

+ Div {2efe>^eWab£,t)ds,.) 
where ds^j, is the local basis of 2-forms on M induced by coordinates. 
Let us define the reduced current £ and the supcrpotential U 

£ = {2ee {fR\ - \fel + Pfey.tbR"\x) e^' - 2V^ {efe'^,) ^';^£,f-^ ds, 

U = 2/' eel^'e:^ ^^^l^^ ds,. = 2e ee^el ^^^^^^^^ ^^^^ 

Accordingly, the Noether current is 

£=£+Div{U) (4.4) 

Notice that the reduced current vanishes on-shell (see equations (3.4)). Conserved quantities 
are then generated by integrating the supcrpotential. 

As usual in gauge-natural theories which are equivalent to a natural theory the correspondence 
among conservation laws is established by means of the so-called Kosmann lift = eJ^e*''V/3^"; 
see [24]. Then the supcrpotential is: 

=2ee (y'^e + Pe^ced^'-'ab-el^^^ bC) ds,, = 

=2ee (V'^e + ^e^e^e«ef e«'''«V^^„) ds^, = (^"^^ 

=2ee^-'e ds^, + 2e/3e''^''«V^Ca ds^, = 

=2eeS7"e ds^, + Div (|e/3e''''^«^„ ds^.p) 
which differs from the Komar supcrpotential computed for the frame e by a pure divergence 
which does not contribute when integrated along closed regions and hence it does not contribute 
to conserved quantities. 

Let us stress that one is forced to choose the Kosmann lift along the new frame e in order to 
obtain a correspondence with the conservation laws of metric-affine formulation. 



(4.3) 
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5. Conclusions and Perspectives 



The result we obtained can be considered from two points of view. Prom the point of view of 
classical gravitational theories, it is interesting to have an equivalent formulation of the usual 
extended theories ,f(R) a la Palatini. Using frames allows coupling to spinors; one has just 
to extend the structure group from S0(?7) to Spin(77). It is also interesting to know that the 
universality theorem extends further to f{^R) theories. Moreover we obtained a non-trivial 
correspondence among conservation laws. This correspondence selects e as preferred frame 
with respect to e. 

From the point of view of Quantum Gravity these models allow a direct approach to quanti- 
zation a la LQG of all f{K) theories; see [25]. Classically these models are known to produce 
modified dynamics for gravitational physics, in particular in Cosmology. Even not considering 
the issue of whether f{R) theories better describe physics than standard GR, it is interesting 
from the theoretical viewpoint to explore the dynamics of such an infinite family of models. This 
could improve the understanding of the emergence of the classical dynamics from the quantum 
world. For example it would be interesting to explicitly see whether the formalism developed in 
these years is able to catch the classical difference of dynamics of these models when compared 
with standard GR. This is particularly relevant in Cosmology where the comparison could help 
in better understanding the relation between LQC and LQG; sec [26]. 

All f{^R) are gauge-natural theories. As such the gauge and diffeomorphism constraint should 
be unchanged. Accordingly nothing should change in defining the Volume and Area operators 
together with their quantizations. Hence the modified dynamics should appear in Hamiltonian 
constraint. It will be interesting, at least from a theoretical viewpoint, to test the proposals for 
quantization techniques against thcrsc extended models; see [27]. 

A forthcoming paper will be devoded to study the Wheeler-DeWitt equation for the extended 
models introduced here. 
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